Universality, maximum radiation and absorption 
in high-energy collisions of black holes with spin 
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We explore the impact of black hole spins on the dynamics of high-energy black hole collisions. 
We report results from numerical simulations with 7-factors up to 2.49 and dimensionless spin 
parameter \ = +0.6, 0, —0.6. We find that the scattering threshold becomes independent of spin 
at large center-of-mass energies, confirming previous conjectures that structure does not matter in 
ultrarelativistic collisions. It has further been argued that in this limit all of the kinetic energy of 
the system may be radiated by fine tuning the impact parameter to threshold. On the contrary, we 
find that only about 60% of the kinetic energy is radiated for 7 = 2.49. By monitoring apparent 
horizons before and after scattering events we show that the "missing energy" is absorbed by the 
individual BHs in the encounter, and moreover the individual BH spins change significantly. We 
support this conclusion with perturbative calculations. An extrapolation of our results to the limit 
7 — y 00 suggests that about half of the center-of-mass energy of the system can be emitted in 
gravitational radiation, while the rest must be converted into rest-mass and spin energy. 
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J. Introduction. Historically, one of the primary moti- 
vations behind the forty-year long effort to develop nu- 
merical relativity was to understand strong-field astro- 
physical phenomena, such as the mergers of black holes 
(BHs) and neutron stars (see e.g. [l|). More recently, 
numerical simulations have begun to shed light on prob- 
lems of fundamental interest in high-energy physics, such 
as trans-Planckian scattering and gauge-gravity dualities 
• A scenario of particular interest in this context is the 
collision of two BHs near the speed of light. This scenario 
has been investigated most extensively in D = 4 space- 
time dimensions for the case of equal-mass, nonspinning 
BH binaries, which are characterized by two parameters: 
the boost factor 7 = (1 — v 2 ) -1 / 2 and the impact parame- 
ter b = L/P, where v is the center-of-mass velocity, L the 
initial orbital angular momentum and P the initial linear 
momentum of a single BH. In the head-on case (b = 0) 
high-energy BH collisions can radiate up to 14 ± 3 % of 
the center-of-mass (CM) energy of the system, and they 
always produce a nonspinning remnant [3j . Grazing col- 
lisions with b =/= 0, on the other hand, result in one of 
the following three outcomes (i) a prompt merger 
for sufficiently small b < b* , (ii) a "delayed" merger for 
b* < b < fescat, or (iii) scattering of the holes to infinity 
for b > & sca t. Here fe sca t denotes the scattering thresh- 
old and b* < b scat the "threshold of immediate merger" : 
by fine-tuning around b* a binary can end up in a near- 
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circular orbit for a time T oc log \ b — b*\, before it either 
separates or merges to form a single Kerr BH Q. In 
[|J, grazing collisions with 7 < 2.9 were found to radiate 
as much as 35 ± 5 % of the CM energy in gravitational 
waves (GWs) and result in near-extremal spins without 
violating the Kerr bound, in agreement with the cosmic 
censorship conjecture (though this may not be the case 
in higher dimensions |6[ ) . A parallel study by Shibata et 
al. [3] investigated the scattering threshold 6 S cat as a func- 
tion of the CM energy. Their simulations of nonspinning 
binaries span boosts up to 7 = 2.3 and yield results that 
are well fitted by 6 scat ~ 2.5 (M/v), where M = 7M0 is 
the CM energy (here and below we use geometrical units 
G = c = 1). Comparisons with BH perturbation theory 
and point-particle collisions in the zero-frequency limit 
provide a satisfactory understanding of the main quali- 
tative features of these simulations [a B 1 but several 
outstanding questions remain. 

Here we perform a systematic analysis of ~ 120 col- 
lisions of spinning and nonspinning BH binaries in four 
dimensions (see e.g. d, [l(| for early results in higher di- 
mensions) to answer two questions of great relevance for 
transplanckian physics and BH production from parti- 
cle collisions, in particular in TeV-scale gravity scenar- 
ios [ll], Ell : i) is the internal structure of the colliding 
objects, here consisting of their spin angular momentum, 
relevant in high-energy collisions? ii) is it possible (as 
suggested in Q) to radiate all of the kinetic energy of 
the system in ultrarelativistic BH encounters? 

Our simulations answer both questions in the negative. 
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FIG. 1. Four waveforms for 7 = 2.49, antialigned spins and 
selected values of b. The b/M = 2.755 case is a triple en- 
counter (two periastron passages followed by a merger). 



We find that spin effects become negligible in the high- 
energy limit, in the sense that both the scattering thresh- 
old and the maximum energy radiated become universal 
functions of 7 (independent of spins). For our largest 
boost (7 = 2.49), grazing encounters never radiate more 
than ~ 60% of the available kinetic energy of the system; 
in fact we find this percentage to decrease for increas- 
ing boost velocity. By monitoring the apparent horizons 
we show that the "missing" kinetic energy is accounted 
for by an increase in the BH mass during the encounter. 
Both observations are important in the context of super- 
Planckian scattering. They justify the use of classical 
general relativity to understand properties of the colli- 
sions and constrain the amount of GWs radiated which 
is a prime example of a relevant property, as this will 
determine the initial mass spectrum of formed BHs. 

Our results thus support from a different viewpoint 
earlier evidence that "matter does not matter" as pro- 
vided in Ref. 



13j . where it was shown that the ultrarel- 



ativistic collision of two bosonic solitons at sufficiently 
high energies leads to BH formation. Similar conclu- 
sions were reached when colliding self-gravitating fluid 
objects [llEl. These results are consistent with hoop 
conjecture arguments, which implies that the nonlinear 
gravitational interaction between the kinetic energy of 
the solitons causes gravitational collapse: the structure 
of the colliding bodies does not matter. 
II. Setup. The simulations presented in this work have 
been performed with the Lean code flal. which is based 
on the Cactus computational toolkit [17| , uses mesh re- 
finement (provided by Carpet [lH) and the apparent 
horizon (AH) finder AHFinderDirect [HHj. Punc- 
ture initial data arc provided by a spectral solver [2l| . 
For obtaining stable evolutions of spinning BHs collid- 
ing with large boosts, we apply two modifications to the 



numerical infrastructure employed in previous studies of 
high-energy BH collisions 0, IS HI: (i) following [H], 
we evolve the conformal factor in terms of the variable 
W = y/x, and (ii) we reduce the Courant factor to 0.45. 

We set up a coordinate system such that the holes 
start on the x axis with radial momentum P x and tan- 
gential momentum P y , separated by a distance d. The 
impact parameter is b = L/P = P y d/P. We extract grav- 
itational radiation by computing the Newman-Penrose 
scalar ^4 at different radii r ex from the center of the col- 
lision. ^4 is decomposed into multipoles ipim using spin- 
weighted spherical harmonics s Yi m of spin- weight s = — 2: 
*4(Ucx,M) = X)j= 2 X)m=-J -2Yi m {9 ,(j>)tpi m (t,r cx ), 
where 9 is measured relative to the x axis. 

Estimates of spurious "junk" radiation inherent in the 
initial data show that it is quite insensitive to the im- 
act parameter, and comparable to our recent findings 
H]; we remove it from reported results in a similar 
manner. Errors due to discretization and finite extrac- 
tion radius are comparable to those reported in [1,|3|- We 
estimate uncertainties in radiated quantities of 3 % and 
15 % for low and high boosts, respectively. The dominant 
error contribution arises from finite resolution, and our 
convergence analysis indicates that our numerical results 
underestimate the radiated energies. 

Contrary to our recent investigation of ultrarclativis- 
tic encounters of spinning, equal-mass BHs in "super- 
kick" configurations [22J, here we expect the dynam- 
ics to be most strongly affected by the "hang-up" ef- 
fect observed for astrophysical BH mergers with spins 
(anti)aligned with the orbital angular momentum j24| . In 
order to investigate the boost dependence of this hang- 
up effect, we evolve three sequences of equal-mass BH 
binaries: (i) a sequence with zero spins, (ii) a sequence 
with spins aligned with the orbital angular momentum 
and dimcnsionlcss spin parameters X = Xi = X2 = 0.6, 
and (iii) a sequence with spins of magnitude \ ~ 0-6 an- 
tialigned with the orbital angular momentum. For each 
sequence we consider four values of the boost parameter 
(7 = 1.22 , 1.42 , 1.88 , 2.49) and for each 7 we simulate 
encounters with about 10 different values of b in order to 
bracket the scattering threshold. 

III. Scattering threshold. As mentioned above, we ex- 
pect a given initial binary configuration to result in either 
a prompt merger, a delayed merger or scatter to infinity. 
Our new simulations confirm this scenario. This is illus- 
trated in Fig. [1] where we plot a subset of representative 
waveforms from the 7 = 2.49 sequence with antialigned 
spins. For small impact parameter (top) the BHs merge 
promptly, and the signal is a clean mcrgcr/ringdown 
waveform leading to formation of a BH with dimension- 
less spin Xf — 0.87. For the second waveform from the 
top the merger is not quite prompt: it shows a pattern 
similar to the scattering waveforms with b/M = 2.772 
shown in the bottom panel, followed by a ringdown. 
The third waveform is an example of delayed merger. 
We chose to display this particular sequence because for 
b/M = 2.755 the binary undergoes a rare triple encounter 
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FIG. 2. Trajectory of one BH from antialigned simulations 
with b/M = 2.755 and 7 = 2.49 (cf. Fig. [TJ. Inset: time evo- 
lution of the irreducible mass Mi rr and of the circumferential 
radius C c of each hole. The circles represent the BH location 
at intervals At = 10 M (corresponding to vertical lines in the 
inset) and have radius equal to the irreducible mass Mi rr . 
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FIG. 3. Critical scattering threshold (upper panel) and max- 
imum radiated energy (lower panel) as a function of 7. Blue 
circles and red crosses refer to the aligned and antialigned 
case, respectively. Black "plus" symbols represent the thresh- 
olds for the four nonspinning configurations studied in this 
paper, complemented (in the upper panel) by results from Q 
for 7 = 1.520. For clarity, we only plot error bars for the 
antialigned-spin sequence. 



consisting of two revolutions (the second close encounter 
is visible as a small "bump" at t/M ~ 550), followed by 
a merger signal with relatively low amplitude. Note that 
the binary radiates and partially absorbs much of the sys- 
tem's kinetic energy during the first scattering encounter, 
therefore subsequent encounters occur at low velocity and 
radiate much less. We graphically display this behaviour 
in Fig. [2j where we plot the trajectory of one BH for the 
configuration b/M = 2.755 and represent snapshots (la- 
beled '1' to '9') of the BH at time intervals At = 10 M by 
circles with radius equal to the irreducible mass. From 
snapshots '2' to '4' we observe a rapid increase in the 
"size" of the black hole; successive snapshots are located 
closer to each other, showing that the BH has slowed 
down during the encounter. 

In order to determine 6 sca t as a function of spin and 
boost we need to distinguish between merging and scat- 
tering collisions. Mergers are easily identified by finding 
a common AH. We identify an encounter as a scattering 
case when the following criteria are met: (i) no common 
AH is found; (ii) the Kretschmann scalar R a p 7 sR a ^ lS at 
the origin approaches zero at late times within numeri- 
cal uncertainties; and (iii) the coordinate trajectories of 
the BHs, Eq. (14) in Ref. [lfl, separate out to values 
comparable to their initial distance. 

The scattering thresholds obtained in this way arc plot- 
ted in the upper panel of Fig. [3l Errors in 6 sca t come 
from numerical truncation error and discrete sampling of 
the parameter space, estimated as follows. In the most 
challenging case (7 = 2.49) our standard-resolution runs 
with grid spacing Ax yield b scat /M = 2.760 for the an- 
tialigned case. By running simulations at two higher res- 
olutions 0.9 Ax and 0.8 Ax, we find b scat /M = 2.741 and 



2.731 respectively, corresponding to about fourth-order 
convergence, a Richardson-extrapolated value b scat /M = 
2.713, and therefore a numerical uncertainty of 0.047. 
The error due to discretization of the parameter space 
(~ 1.6 x 10~ 3 ) is negligible in comparison, so we adopt 
8b sca t ~ 0.05 as a conservative error estimate. Fig. [3] 
shows that the scattering threshold is spin-independent 
in the limit 7 — > 00. Our nonspinning simulations are 
consistent with the results obtained by Shibata et al. Q 
at lower boosts. They are also compatible with the shock- 
wave analysis in Table II of (25|, that suggests a critical 
impact parameter b sca t/M > 1.68 in the ultrarelativistic 
limit. Note that the shock-wave construction overesti- 
mates the total radiated energy in head-on collisions by 
a factor ~ 2 Q; if this is also the case for grazing col- 
lisions, one would expect that shock-wave methods may 
underestimate b sca t/M by a similar factor. 
IV. Maximum radiation. As pointed out in [H, [26| , 
the total energy £" ra( j/M radiated in grazing BH col- 
lisions increases steeply as the impact parameter ap- 
proaches b* or 6 sca t. This is true also for spinning bi- 
naries. For 7 = (1.22, 1.42, 1.88, 2.49), respectively, we 
find the maxima in E la d/M plotted in the lower panel 
of Fig. [3J For reference, the initial fraction of total en- 
ergy in the form of kinetic energy K/M — (7 — l)/7 
is (17.9,29.4,46.8,59.8)%, and for the spinning cases we 
have (4.2, 3.6, 2.6, 1.7)% in initial spin energy. 

For a subset of scattering runs where we monitored the 
apparent horizon as a function of time, we list estimates 
of radiated energies and spins before/after scattering in 
Table HI The accuracy of the radiated energy estimates 
is limited by discretization errors, the presence of spuri- 
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ous initial radiation and by our discrete sampling of the 
impact parameter space, but the numbers clearly reveal 
two striking features: (i) the maximum radiated energy 
varies only mildly with spin at any given 7, and (ii) for 
mild boosts the maximum energy radiated is comparable 
to the initial kinetic energy of the spacetimes; however as 
7 increase the ratio drops, down to ~ 60% for 7 = 2.49. 
This observation prompts two questions. Where has the 
remaining kinetic energy gone? Why does the deficit in- 
crease with boost? 
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TABLE I. Fractional kinetic energy radiated and absorbed 
for representative numerical simulations with aligned spins 
(f), antialigned spins (\) or nonrotating BHs ('0'). See the 
main text for definitions of the various quantities. We also 
list spin estimates Xi and Xs before and after scattering. Xi 
is measured at a time ~ 20 M after the beginning of the sim- 
ulation, and small deviations from the initial data parameter 
j — ±0.6, can presumably be attributed to the BHs absorb- 
ing an increasing amount of junk radiation as 7 increases. 

V. Absorption. The answer to these questions is found 
in the apparent horizon dynamics of the individual holes 
before and after the first encounter. We have analyzed 
the data in detail for a set of scattering events where 
the individual holes separate sufficiently after first en- 
counter to warrant application of the isolated horizon 
limit. Specifically, we measure the equatorial circum- 
ference C = A-nM and the Christodoulou mass M? T of 
each BH before and after the encounter. The inset of 
Fig. [2] shows the variation of these quantities with time 
in a typical simulation: absorption typically occurs over 
a short timescale ~ 10M. Since the apparent horizon 
area ^ AH = WttM^ = [C c 2 /(2tt)](1 + ^Jl-x 2 ), in this 
way we can estimate the rest mass and spin of each hole 
before (Mi, x0 an d after (M s , Xs) the first encounter. 
We define the absorbed energy i5 abs = 2(M S - Mi). The 
results in Table |T] show that the sum (E ra< i + E a b s )/M 
accounts for most of the total available kinetic energy 
in the system, and therefore the system is no longer 
kinetic-energy dominated after the encounter. A fit of 
the data yields E IaA /K = 0.46(1 + IT/7 2 ) and E ahs /K = 



0.55(1 — I/7), suggesting that radiation and absorption 
contribute about equally in the ultrarelativistic limit, and 
therefore that absorption sets an upper bound on the 
maximum energy that can be radiated. 

The fact that absorption and emission are comparable 
in the ultrarelativistic limit is supported by point-particle 
calculations in BH perturbation theory. For example, 
Misner et al. [27[ studied the radiation from ultrarela- 
tivistic particles in circular orbits near the Schwarzschild 
light ring, i.e. at r = 3M(1 + e). Using a scalar- 
field model they found that 50% of the radiation is ab- 
sorbed and 50% is radiated as e — > 0. We verified by 
an explicit calculation ignoring self-force effects that the 
same conclusion applies to gravitational perturbations of 
Schwarzschild BHs (cf. [28j). A recent analysis including 
self-force effects finds that 42% of the energy should be 
absorbed by nonrotating BHs as e — > (cf. Fig. 4 of |{|). 

Rather than considering particles near the light ring, 
we can model our problem using particles plunging ultra- 
rclativistically into (for simplicity) a Schwarzschild BH. 
Davis et al. [29j first computed the energy absorbed when 
a particle of mass m falls from rest into a Schwarzschild 
BH of mass Mbh, and they found the remarkable result 
that the total absorbed energy (summed over all multi- 
poles £ > 2) diverges. Physically, the divergence is due to 
the fact that most of the absorption occurs near the hori- 
zon, so we must go beyond the point-particle approxima- 
tion and introduce a physical cutoff at £ max ~ 7tMbh/2ot 
to take into account the finite size of the infalling par- 
ticle. For comparable-mass encounters it is reasonable 
to truncate the sum at £ = 2. By adapting the BH 
perturbation theory code of [3(|, we extended the cal- 
culation of (29| to arbitrary spacetime dimensionality 
and to generic particle energies po = E/m. Details of 
this study will be reported elsewhere. For D = 4, and 
truncating the sum at I = 2 as in [2!| when we com- 
pute absorption, our calculation shows that the radiated 
(absorbed) energy is i^d.abs = fc rad,abs(pom) 2 /M B n, 
with fc rad = (1.04 x 10" 2 ,3.52 x 10~ 2 , 0.119, 0.262) and 
fc abs = (0.304, 0.310, 0.384, 0.445) for Po = (1, 1.5, 3, 100), 
respectively. So, again, in the UR limit the point-particle 
model predicts a roughly comparable amount of emis- 
sion and absorption. As a consequence of this significant 
energy absorption, in the large-7 limit close scattering 
encounters between two arbitrarily small (in rest mass) 
BHs can result in two, slowly moving BHs with rest mass 
increased by a factor of order 7. Figure [2] provides an ex- 
plicit example of this scenario, where BHs grow in size 
due to absorption and slow down after scattering. 

Another remarkable implication of Table U is that close 
high-energy scattering encounters can drastically modify 
the spin magnitude, presumably due to some combina- 
tion of absorption and dynamical tidal interactions. For 
example, when 7 = 2.5 the BH spins decrease from ~ 0.6 
to ~ 0.3 in the aligned case, from ~ 0.6 to ~ 0.1 in 
the antialigned case, and we measure a post-scattering 
spin x ~ 0-2 for initially nonspinning encounters. These 
changes in dimensionlcss spin parameter correspond to 
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roughly the same total angular momentum being trans- 
ferred to the black holes during the interaction, indepen- 
dent of the initial spin. 

VI. Ultrarelativistic extrapolation and a conjec- 
ture. We can now try to estimate the maximum radiated 
energy when 7 — > 00 as follows. In Ref. Q, we estimated 
that for head-on collisions (b = 0) we can radiate at most 
E xa d/M ~ 0.14 in this limit. For each Lorentz factor 
7, the increase in radiation induced by fine-tuning near 
threshold can be characterized by the ratio 72 = E Ya d(b = 
0)/i?rad X - This ratio increases with boost: for our largest 
7 it is approximately 72 ~ 0.2, and a fit to our data yields 
72.(7) = 0.34(1 — I/7). Using this fit in combination with 
results from Ref. [| we find E^/M w 0.14/0.34 - 0.41 
as 7 — > 00. A more conservative upper bound on the 
total radiated energy at threshold can be obtained as- 
suming that 72 is a monotonically increasing function of 
7, but using the last data point in our simulations as a 
lower limit on 72: this yields a very conservative upper 
limit E^/M ^ 0.14/0.2 = 0.7 as 7 -> 00. These results 
are consistent (within the errors) with our discussion of 
Table HI which indicates that radiation in high-energy en- 
counters accounts for roughly 0.46 of the available energy, 
the rest being absorbed. Therefore our simulations seem 
to settle the long-standing question of whether it is possi- 
ble to release all of the CM energy as GWs in high-energy 
BH collisions: the answer is no. 

Two crucial assumptions underlie the study of BH pro- 
duction from high-energy particle collisions. The first as- 
sumption, that BHs are indeed produced as a result of 
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